o 
o 

(N 

-t— > 
o 

O 



o 
o 

u 

Oh' 



X3 
O 

o 



> 
m 

(N 
O 



o 



o 



X 
S3 



Anisotropy of Vortex-Liquid and Vortex-Solid Phases in Single Crystals of 
Bi2Sr2CaCu208+5: Violation of the Scaling Law 

J. MirkovicQ, S. E. Savel'ev, E. Sugahara, and K. Kadowaki 
Institute of Materials Science, University of Tsukuba, Tsukuba 305-8573, Japan 

The vortex-liquid and vortex-solid phases in single crystals of Bi2Sr2CaCu2 08+i placed in tilted 
magnetic fields are studied by in-plane resistivity measurements using the Corbino geometry to 
avoid spurious surface barrier effects, ft was found that the anisotropy of the vortex-solid phase 
increases with temperature and exhibits a maximum at T « 0.97 Tc. In contrast, the anisotropy of 
the vortex-liquid rises monotonically across the whole measured temperature range. The observed 
behavior is discussed in the context of dimensional crossover and thermal fluctuations of vortices in 
the strongly layered system. 



The anisotropic static and dynamical properties of su- 
perconductors have commonly been described by the 
three-dimensional Ginzburg-Landau (hereafter 3DGL) 
theory with a single parameter of anisotropy 7 defined 
as a ratio {rric/mab)^^'^ of the eflFective masses TOc and 
niab along the c-axis and in the a6-plane, respectively [Q. 
The anisotropy is related to the basic parameters of a 
superconductor as 7 = £,ah/ic = K/Kb = Hc2\\/Hc2±, 
where ^ab and are the in-plane and the out-of-plane 
coherence lengths, Ac and Xab are the magnetic field pen- 
etration depths along the c-axis and in the a6-plane, while 
Hc2± and i?c2|| are the out-of-plane and the in-plane up- 
per critical magnetic fields. 

Based on 3DGL, Blatter et al. |^ obtained the general 
scaling law, which can be applied for the description of 
different physical properties of anisotropic superconduc- 
tors in oblique magnetic fields. For instance, the resis- 
tivity in a magnetic field H , tilted away from the c-axis 
for the angle 6, depends only on the reduced field as 



p{H, 9) = p{H J cos^ 9 + sin^ 6I/72). 



(1) 



The scaling law for the vortex lattice melting magnetic 
field was derived in a similar manner as 



iJ™^"(6') = iJ-™<='*(6l = 0)/-y/cos2 9 + sin^ 9h^. (2) 

In previous studies it was found that the scaling 
law describes well experimental data such as resistivity 
[^,||, magnetization and thermodynamic measure- 
ments in YBa2Cu307_5 with anisotropy parameter 
7 « 5 — 7 [0 . Concerning the higher anisotropic layered 
systems, it was recognized earlier that the dissipation in 
Bi2Sr2CaCu208+(5 depends practically only on the c-axis 
magnetic field component except very close to the 
a6-plane which is in agreement with the 3DGL scal- 
ing law with high anisotropy parameter. However, it was 
reported [lO| ] that the resistivity data near the a5-plane 
agree better with Thinkam's thin film (hereafter 2DT) 
model jl^], than with the 3DGL theory. Later, it was 
also found that the measured anisotropy parameter for 
Bi2Sr2CaCu208+5 thin films increases with temperature, 
according to the 2DT model, followed by an indication of 
the two to three-dimensional (2D-3D) crossover near the 



transition temperature [|l^,|l^ . On the other hand, re- 
cent experimental studies of the first-order vortex-lattice 
melting transition (hereafter VLMT) have shown a linear 
decay of the out-of-plane melting field with an increasing 
in-plane magnetic field, even near the c-axis in a 

strong contrast to the above mentioned approaches. In 
a strict sense, it is reasonable to expect that the scal- 
ing law based on the 3DGL model may not be appli- 
cable to the extremely anisotropic layered system such 
as Bi2Sr2CaCu208+5. In this material, the coherence 
length along the c-axis at low temperatures is smaller 
than the lattice parameters of the crystal, which means 
that the layerness of the material has to be considered 
Thus, the 3DGL scaling approach based on a 
continuous medium, is not valid due to the discreteness 
of the system, and therefore, it is not at all trivial that 
the expressions (|l|) and (||) are still valid. 

Hence, there is still a question of how the anisotropy 
manifests itself in the different vortex phases in 
Bi2Sr2CaCu208+5 single crystals. In this Letter, we 
report the comparative study of the anisotropy of the 
vortex-solid and vortex-liquid phases in single crystals of 
Bi2Sr2CaCu208+5. The anisotropy of the first-order vor- 
tex lattice melting transition increases with temperature 
in the range of T < 0.97 Tc, whereas at higher tempera- 
tures, the VLMT anisotropy decreases. In contrast to the 
anisotropy of the VLMT, the anisotropy of the in-plane 
resistivity in the vortex liquid phase increases monoton- 
ically with temperature in the whole measured range. 

The in-plane resistivity measurements were performed 
for three as-grown single crystals of Bi2Sr2CaCu208-i-5 



with transition temperatures T 



{I) 



90.3 K, T, 



11) 



86.0 K, and ri^^^^ = 90.0 K, for the samples #sl, #s2 
and #s3, respectively. In order to probe the true bulk 
properties, i.e., to avoid the surface barrier effects which 
occur in conventional strip geometry [p6[|, we have used 
the Corbino electrical contact configuration (see inset 
in Fig. la). The diameters of the Corbino discs were 
i:>(^) = 1.9 mm, L>(^^) = 1.95 mm, and D^"^'^ = 2.7 
mm while the thickness was t k, 2Q fim for all samples. 
The resistance has been measured by using the ac lock- 
in technique at a low frequency of 37 Hz. The magnetic 
field, generated up to 70 kOe by a superconducting split 
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FIG. 1. (a) The magnetic field dependence of tlie re- 
sistance at various field orientations from H || c (0°) to 
H II ab(90°) at r=86.6 K (sample #sl). Inset in (a): electri- 
cal contacts in the Corbino geometry, (b) The — Hat phase 
diagram of the vortex lattice melting transition at T = 86.6 
K: the filled symbols mark the first-order phase transition, 
while the open ones correspond to the continuous resistivity 
transition. Inset in (b): the part of the VLMT phase dia- 
gram near the c-axis, plotted at a magnified scale, (c) The 
equi-resistance contour at 86.6 K (open symbols) in the vor- 
tex-liquid phase. The melting transition from (b) is replotted 
at the centre. The solid line is the ellipse with anisotropy 
value of 7l ~ 155, while the dotted line corresponds to the 
phenomenological formula wih 'yphen = 200. 



magnet, was rotated by a fine goniometer with angular 
resolution of 0.01°. 

A typical set of the resistance curves is presented in 
Fig. la for T = 86.6 K. The sharp resistance drop, at- 
tributed to the VLMT ||l^ , is clearly detected across the 
wide angular range {9 < 89.86°) ||l5|. Using the resis- 
tance criteria of i? = 10~^ri, which is set somewhat below 



the resistivity kink, and R = 10~^fl, which is well above 
the VLMT anomaly, we construct the phase diagram of 
the VLMT and the equi-resistance contour for the vortex- 
liquid in the He — Hab phase plane (Fig. l(b,c)). The 
VLMT phase Hne H^'=^\Hab) (Fig. lb) exhibits the pe- 
culiar step- wise shape Namely, starting from the 
c-axis, the out-of-plane component of the VLMT, 
decays linearly with increasing Hab p4[|, w hich is associ- 
ated with the crossing vortex latticep^ . As the mag- 
netic field further approaches the a6-plane, the linear de- 
pendence H™'^^^{Hab) sharply transforms to the plateau 
(inset in Fig.l b), which is attributed to the trapping of 
pancake vortices by Josephson vortices Very close 
to the a6-plane {6 « 89.96°), the first-order phase tran- 
sition (filled symbols in Fig. lb) changes to the continu- 
ous resistance transition (open symbols), which exhibits 
the cusp in the He — Hab phase plane In strong 
contrast to the linear decay of H™°^^{Hab) around the 
c-axis, the equi-resistance line Hf {Hab) (Fig. Ic) in the 
vortex-liquid phase demonstrates smooth quadratic be- 
havior. Near the a6-plane, the equi-resistance contour 
still has a tail, as if it traces the VLMT cusp. The cusp- 
like shape of both the VLMT line and the equi-resistance 
contour, pronounced near the a6-plane, is possibly asso- 
ciated with the intrinsic pinning which, surprisingly, 
seems to be active even in the vortex liquid phase, quite 
far from the VLMT. 

The next question is how to estimate the anisotropy 
of the vortex liquid phase and the VLMT. The smooth 
quadratic dependence H^{Hab), can be derived from 
the general symmetry law p{Hc, Hab) = piHc,-Hab), 
assuming that the resistivity is an analytical function 
of Hab in the vortex-liquid. The simple analysis pl| ] 
gives the angular dependence H^{Hab) = H^{0) — 
^ab/ i^c i^)lphen) ^O'" limited angular range 0° < 
9 < 90° — 57.3° /'jphen with phenomenological anisotropy 
parameter ^pheniT, He). This equation succesfuUy fits 
our data at 86.6 K in the vortex liquid phase except 
very close to the G6-plane with ^phen ^ 

200 (dotted Hne 

in Fig. Ic). It is worth noting that the phenomeno- 
logical angular dependence H^{Hab) coincides with the 
angular dependence of the upper critical field Hc2{9) 
for thin superconducting films obtained by Tinkham 
as H,2[9)cos9/H,2±+Hl^{9)sin^9/-fl^^Hl^^ = 1 
with the field independent anisotropy parameter 72DT oc 
1/^1 — T /Tc- On the other hand, to describe the equi- 
resistance contours in the vortex liquid phase, it is conve- 
nient to apply the more common 3DGL model. Following 
the 3DGL theory, the H^{Hab) line should be an ellipse 
{H^{Hab)f + HlJ^l^aL = {H^f with anisotropy 
l^DGL being independent of magnetic field and temper- 
ature. Fitting the curve in Fig. Ic to an ellipse with an 
anisotropy parameter of 155, we also find a reasonable 
agreement only in a limited angular range 9 < 89.75°. 
Moreover, the empirical anisotropy 7/, of the vortex 
liquid, obtained by the elliptical fitting of the equi- 
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FIG. 2. The temperature dependence of 7l, extracted 
by fitting tiie equi-resistance lines to an ellipse, for samples 
T^sl and #s2. The dashed line corresponds to the depen- 
dence 7L = 28/(1 - T/Tcf"^. Inset: The dotted curves 
show the R{H) data obtained at r=89.5 K at various mag- 
netic field orientations with respect to the c-a:xis (from left, 
6 = 66°, 74° , 80° , 84° , 86.5° , 88.05° , 88.75° , 89.2° , 89.39° , 89.6° ) . 
The solid curves represent scaled resistance curves according 
to equation (j^) with 7 = 71, = 350. 

resistance contours, increases strongly with temperature 
in the whole measured temperature interval (Fig. 2), 
which indicates that the 3DGL scaling law is violated. 
Thus, the anisotropy parameter "/l extracted from the 
ellipse can no longer be attributed to the ratio of the ef- 
fective masses as assumed in the 3DGL theory and may 
reflect the two-dimensional character of the vortex sys- 
tem in analogy to the 2DT approach jllj. Interestingly, 
the resistivity curves above the VLMT kink, measured 
at different field orientations {9 < 89.6°, T = 89.5 K, 
R < 2 ■ collapse into a unique one if the curves 

R{H, 9) are scaled in a common way by using eq. (|l]) (see 
inset in Fig. 2) with 7 = 71,. Hence, the anisotropy 
depends weakly on the resistance level in the vortex liq- 
uid but changes sharply around the vortex-lattice melt- 
ing transition (see Inset in Fig. 2), implying different 
dimensionality of the probed vortex phases. 

Recognizing the fact that the VLMT does not follow 
the 3D scaling law (H), the effective anisotropy of the 
VLMT can be estimated by the ratio 7me;t = H*fj/H* 
(for definition of H* and i?*^, see Fig. lb). It is neces- 
sary to emphasize that the VLMT phase lines, obtained 
at different temperatures, fall roughly on a unique line if 
plotted in the phase plane Hc/H* — Hab/H*^, as shown 
in Fig. 3a. Therefore, the whole VLMT line in the 
He — Hah plane changes proportionally with temperature, 
which justifies the chosen definition of the anistropy of 
the VLMT. The temperature dependences of the 
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FIG. 3. (a) The phase diagrams of the vortex melting tran- 
sition plotted in the phase plane Hc/H* — Hab/H*i, at differ- 
ent temperatures for #sl. The left and right insets display 
the temperature dependence of and H* , respectively (for 
definition of the fields and H*, see Fig. lb), (b) The 
temperature dependence of the VLMT anisotropy defined as 
77neit = H*fj/ HI, obtained for three samples. The solid line is 
the empirical equation for 'ymeitiT) discussed in the text. In- 
set: the calculated temperature dependence of the anisotropy 
in the frame of the model [23] . 

characteristic fields * and iJ*^ exhibit different curva- 
tures (see insets in Fig. 3a), which result in a maximum 
of the temperature dependence of the VLMT anisotropy 
ImeitiT) at T w 0.97 as it is seen in Fig 3b. The 
obtained temperature dependence of ^meit may be ap- 
proximated by the empirical equation ^meit{T) — 70 (1 ~ 
7/Tci)^/(l - T/Tc2)" (the solid line in Fig. 3b) with 
parameters = 89.4 K, 7^2 = 90.7 K, a = 1.08 
and P = 0.6 obtained by fitting the fields H*{T) = 
H*{0){1 - T/T^^r and Hl,{T) = H*jO)il - T/T^.f 
for the sample #sl. The two characteristic tempera- 
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tures (Tci, Tc2) in the above empirical formula, may 
appear due to the fact that the thermodynamical vor- 
tex fluctuations have different influence on the interlayer 
and the intralayer-superconductivity in the tilted mag- 
netic fields. The similar maximum of 7(T), accompa- 
nied by the disappearance of the cusp in the dependence 
Honset{0) {R{Honset) = 0) near the a5-plane, was ob- 
served earlier in thin films of Bi2Sr2CaCu208+5 |Q, al- 
though there was neither indication of the VLMT nor 
clear evidence which vortex phase was probed. The phe- 
nomenon was interpreted ||T^ as a crossover from 
2D to 3D behavior, which occurs when the coherence 
length exceeds the distance between Cu02 layers. In 
agreement with this scenario, the cusp in the He — Hah 
phase diagram, associated with the continuous resistiv- 
ity transition, also disappears around the temperature 
(see Fig. 3a) at which jmeitiT) is at a maximum [ p2[ . 
However, based on the 2D-3D crossover, it is impossi- 
ble to explain why the anisotropy of the vortex liquid 
phase still increases at temperatures T > 0.97 Tc where 
lmeit{T) decreases. Moreover, the 2D-3D scenario pro- 
posed in [Tsl describes only the temperature and angular 
dependence of the upper critical field Hc2iO,T), which, 
strictly speaking, can not be applied to the VLMT and 
the equi-resistance contours. 

Another possible origin of the temperature depen- 
dence of the anisotropy could be related to the ther- 
mal fluctuations of vortices, which suppress the Joseph- 
son coupling in-between the Cu02 layers more efficiently 
than the superconductivity within the layers. As a 
result, Ac(T) increases faster than XabiT) and, thus, 
the anisotropy 7 = Xc/Xab rises with increasing tem- 
perature. According to the model proposed for the 
magnetic field parallel to the c-axis [Q, the temper- 
ature dependence of the anisotropy is determined by 
the equation 7^/79 = exp (S{Hc,T)^^ /^q) jH], where 
6 = TTTHcsXlf,{T)-/o/2<i>l with 70 = 7(T = 0). In order 
to take into account the temperature dependences of the 
VLMT and equi-resistance magnetic fields, the empirical 
equation = Hq{1 — T/Tc)°' is used. Since 7 is a mono- 
tonical function of (5 oc T(l — T/Tc)°'~^ , the tempera- 
ture dependence of 7 has a maximum approaching Tc, if 
a > 1, otherwise 7 monotonically increases with T. The 
inset in Fig. 3b displays the temperature dependence of 7 
obtained by using the above equations with 70 = 300 and 
a = 1.08, corresponding to the measured VLMT (inset 
in Fig. 3b), and a = 1, as assumed for the vortex- liquid 
phase. Nevertheless, the calculated value of the magnetic 
field Hq — 29 kOe disagrees with the experimental value 
of the c-axis field component, but is close to the total 
applied magnetic field. The discrepancy could be related 
to the presence of the Josephson vortices, which have 
not been considered in the model . The pancake vor- 
tex sublattice interacts with the Josephson vortex sub- 
lattice in the crossing lattice structure (set in the tilted 



magnetic fields) providing the additional mechanism of 
the renormalization of the anisotropy in the vortex-solid 
phase [|8|jl|. 

In summary, in contrast to the prevailing belief, we 
have found clear experimental evidence that neither the 
3DGL scaling law nor the 2D Tinkham's model consis- 
tently describe the resistivity in the vortex-liquid phase 
as well as the behavior of the vortex lattice melting tran- 
sition, which can qualitatively be accounted for by the 
thermal fluctuations of vortices. 

We thank M. Tachiki, L. N. Bulaevskn, A. E. Koshelev, 
and V. M. Vinokur for stimulating discussions. 
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